Harmonic entanglement with second-order non-linearity 
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We investigate the second-order non-linear interaction as a means to generate entanglement between fields 
of differing wavelengths. And show that perfect entanglement can, in principle, be produced between the 
fundamental and second harmonic fields in these processes. Neither pure second harmonic generation, nor 
parametric oscillation optimally produce entanglement, such optimal entanglement is rather produced by an 
intermediate process. An experimental demonstration of these predictions should be imminently feasible. 
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Second-order nonlinear processes have found many appli- 
cations in electronics, mechanics, optics, as well as numer- 
ous other fields. More recently, in the field of quantum op- 
tics it has been discovered that they offer an ideal mechanism 
for the generation of entangled states of light. As a result 
many fundamental theories of quantum mechanics have been 
tested at an unprecedented level 1 1]. It is surprising then to 
find that, to date, a comprehensive treatment of entanglement 
generated through second-order non-linear processes has not 
been presented. Here, we perform a broad analysis of entan- 
glement in degenerate processes of this kind. We demonstrate 
that, in principle, perfect entanglement can be achieved be- 
tween the fundamental and second harmonic fields involved 
in the process. This entanglement between harmonically re- 
lated fields is referred to as harmonic entanglement. It is well 
known that quantum correlations and harmonic entanglement 
can be achieved via second harmonic generation (SHG) Q], 
and furthermore that highly squeezed states of light can be 
obtained in the complimentary process, optical parametric os- 
cillation (OPO) 010]. However, we find that optimal har- 
monic entanglement is not generated by either of these pro- 
cesses, but rather by an intermediate process, pump depleted 
optical parametric amplification (OPA). In contrast, previous 
investigations of the quantum optical properties of OPA have 
almost exclusively been limited to the small region of param- 
eter space where pump depletion is insignificant. 

The study of entangled states began as a means to test the 
counter intuitive predictions of the theory of quantum me- 
chanics. In recent years, however, focus has shifted as a result 
of the realization that non-classical states of light can enhance 
many measurement, computation, and communication tasks 
0,1^. That quantum mechanics, in principle, allows such en- 
hancement is significant in and of itself, and has motivated a 
great deal of theoretical work. However, emphasis has also 
been placed on experimental demonstrations, which as yet 
have been primitive by comparison. New tools are required 
for significant experimental progress to be made. Both SHG 
and OPO have found many applications in quantum optics ex- 
periments 0,1^ 01, becoming standard tools in the field. The 
pump depleted OPA studied here has the potential to become 
equally significant. The generation of strong entanglement 
between optical fields at vastly different frequencies would 
immediately facilitate many interspecies quantum information 




FIG. 1: Schematic for the generation of optical harmonic entangle- 
ment. 



protocols, for example the frequency of states of light could 
be drastically altered using interspecies quantum teleportation 
protocols. Such a protocol would enhance the integrability 
of disparate nodes in a quantum information network. Ap- 
plications could be envisaged in any situation where a non- 
classical link is required between two experiments at differing 
optical frequencies. Such a situation might arise for example 
in experiments involving two different atomic species, or if a 
connection is required to an atomic frequency standard. 

Analysis: The system under analysis in this paper consists 
of a second order nonlinear medium enclosed within an opti- 
cal resonator as shown in Fig.^ The resonator is coupled to 
the environment through two partially reflective mirrors. One 
mirror represents an input/output coupler, while the other rep- 
resents uncontrollable coupling loss. The non-linear medium 
induces an interaction between the two intra-cavity fields, giv- 
ing LiOJ 

a = -K.aa + ea% + A\n 
h — —Kbb—^ea'^ + Bin, 

where a and b are Heisenberg picture annihilation operators 
describing the intra-cavity fundamental and second harmonic 
fields respectively; Ka and Kb are the associated total res- 
onator decay rates; e is the nonlinear coupling strength be- 
tween the fields; and Ain and Bin represent the accumulated 
input fields to the system. Throughout this paper the partially 
reflective mirrors modeling input/output coupling and loss are 
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FIG. 2: A plot of the EPR entanglement strength as a function of the 
fundamental, a, and second harmonic, (3, driving field amplitudes. 
Bi-stable OPA (i) above and (ii) below threshold. Parametric de- 
amplification regime (iii) above and (iv) below threshold. 



distinguished with the subscripts '1' and '2', respectively; 
while the input and output fields are denoted by the subscripts 
'in' and 'out'. Using this terminology — Kai + Ka2^ 

Kb ~ Khl + Kh2, A;n = 



V2Kal^l,in + V2'«a2^2,in, and 

\/2Khii?i.iii + \/2Kb2_B2,in- The solution to Eqs. ([Q 
is obtained through the technique of linearization, where op- 
erators are expanded in terms of their coherent amplitude and 
quantum noise operator, so that a mode a; — a + Sai with 
{5ai) = 0, and second order terms in the quantum noise oper- 
ators are neglected. First, the steady state coherent amplitudes 
of the fundamental and second harmonic intra-cavity fields are 
obtained. After a stability analysis has been made on the solu- 
tions, the system shows a range of interesting behavior such as 
mono-stability, bi-stability, out-of-phase mono-stability, and 
self-pulsation 1 10|. For clarity, we normalize the driving fields 
to their respective critical transition values, i.e. the critical am- 
plitude, ai,in,c, for self-pulsation in SHG, and the threshold 
amplitude /3i in,c for OPO, such that 



{2Ka + Kb)[2Kb{Ka + Kb)]^^^ 



/3l,in,c 



e^/2Kbi 



(2) 
(3) 



The normalized driving fields are then denoted = 
aian/aian,c and /3d = /3i,in//3ian,c- The quantum fluctua- 
tions of the intra-cavity fields can be obtained from the fluctu- 
ating part of Eqs. Q 



5h = —Kb6b — eaSa + Bm 



(4) 



These equations can be easily solved by taking the Fourier 
transform into the frequency domain. We want to investigate 
harmonic entanglement between the amplitude X+ and phase 



X~ quadratures of the two output fields. These quadratures 
are related to the annihilation operators via = a + and 
X^ — i{a'^ — a). Writing the solution to Eqs. (|4} in terms of 
field quadratures, we find 
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where {X^,Xj^} and {X~^ j^^j- are me intra-cav 

accumulated input field quadratures, respectively; and 



Xg in} ^1"^ '■^^ intra-cavity and 
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Ka-iuj + e|/ 
— e|/3| sin 
— ejal cos 9a 
— e|a| sin 0^ 

Kb — ILO^ 



I cos Ufj 
I cos 0f3 



(6) 
(7) 
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with 9a = Arg(a), Op = Arg(/3), and uo is the side-band de- 
tection frequency. The fields output from the resonator can 
then be directly obtained using the input-output formalism, 

xi 



U ,out - V2KQl^^-^Aian' ^Bl.out - ^'^'^blX^-^^-Q^ -^^ 

11 ill . We can now proceed to investigate the presence of en- 
tanglement between the output fundamental and second har- 
monic fields. A bi-partite Gaussian entangled state is com- 
pletely described by its correlation matrix lll2ll which has the 
following elements 



X? 



Chi 



1 



mn 



+ X,',X,^„)-(X^)(Xi) (12) 



where {k,l} e {+,-} and {m,n} g {Ai^out, Bi,out}- The 
standard characterization of continuous variable entanglement 
is to measure the quantum correlations between two fields via 
the EPR criterion fl3ll and to apply the inseparability crite- 
rion II12I1 . Before the inseparability criterion can be applied 
however, the correlation matrix is required to be in standard 
form II, which can be achieved by application of the appropri- 
ate local-linear-unitary-Bogoliubov-operations (local rotation 
and squeezing operations) |12]. The product form of the de- 
gree of inseparability fTil is given by 



CfCj 



fc2 -f l/P 



where 



Cf = . 



k = 




±± 

bb 



±± 



1 



(13) 



(14) 



(15) 



X < 1 is a necessary and sufficient condition of inseparability 
and therefore entanglement. We will use X as a measure of en- 
tanglement. The degree of EPR paradox, which measures the 
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FIG. 3: Harmonic entanglement observed in (a) SHG and (b) OPO. 
Botii tiie EPR (solid lines) and the inseparability (dotted lines) en- 
tanglement measures are plotted as function of the fundamental and 
second harmonic field amplitudes, respectively. 



level of apparent violation of the Heisenberg uncertainty prin- 
ciple achieved by the state, can also be found from elements 
of the correlation matrix. 



-f 



r'++ 
^bb 



n — \'^ab P 



C, 



66 



(16) 



The state is entangled when e < 1. Note that e is minimized 
when the correlation matrix is in standard form. We there- 
fore restrict our analysis to this case, denoting the optimized 
degree of EPR paradox as Zo- In this work, the optimization 
of the correlation matrix was achieved numerically. The pa- 
rameters used for making our calculations were = 1 , 
Ka2 — 0.01 , Kbi = 10 , Khi = 0.1 , e = 1 , w = which 
were chosen to resemble the squeezing/entanglement genera- 
tors of recent experiments . 

Results: Using these analytical results, we mapped out the 
degree of optimized EPR paradox (Eq) of harmonic entangle- 
ment across the parameter space of fundamental and second- 
harmonic driving field amplitudes. This has been visualized 
in Fig. 12] where a greyscale-compatible color map has been 
assigned to a range of EPR values on a dB log (base 10) scale. 
Areas that are printed in dark ink, signify strong entangle- 
ment. White areas are not entangled, and are either at the 
standard quantum limit or greater (classical thermal states). 
The vertical {(3d) and horizontal (ad) axes of the plot coin- 
cide with the special cases of OPO and SHG respectively. Re- 
gions of interest have been individually labeled. Note that in 
order to represent bi-stability in the map, the left and right- 
hand sides of the map must be mentally folded along the 
OPO axis. The map of driving fields is clearly divided into 
the three distinct areas that correspond to the stable branches 
of the steady-state solutions. These are; (i) the bi-stable re- 
gion above OPO threshold, (ii-iii) the mono-stable solution 
showing parametric amplification at (ii) and de-amplification 
at (iii). And finally, the so-called out-of-phase solution, which 
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FIG. 4: Polar plot of harmonic entanglement strength with a total in- 
put power of (a) 90% of threshold and (b) 400% of threshold. Shaded 
areas denote the presence of harmonic entanglement. 



has a complex-valued fundamental field solution, is labeled 
with (iv). We find that the system produces a maximum en- 
tanglement of up to Eo = 13.8 dB at the intersection of (i-ii) 
and also when crossing (iii-iv). These regions correspond ex- 
actly to the critical boundaries where a branch change in the 
classical solution occurs. 

A more detailed analysis of entanglement in SHG is shown 
in Fig. |3ja) where the EPR measure is plotted concomitant 
with the inseparability criterion as the amplitude of the driv- 
ing field is varied. The strength of entanglement finds a max- 
imum of So = 1.9 dB and T ~ 3.0 dB for the system driven 
with a fundamental field of ad — ±0.1. In SHG, entangle- 
ment is produced for all non-zero driving field amplitudes. 
However, the strength of the entanglement does not become 
arbitrarily high if loss in the system is made arbitrarily small. 
This is a bound set by the SHG process itself. Note that for a 
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traveling-wave SHG process the strength of entanglement has 
been found to be Eq ~ 12 dB t2J- Contrast this behavior with 
the OPO process, which for below threshold, is not entangled, 
and only becomes entangled once the system is pushed be- 
yond threshold. This is shown in Fig. |3jb). The maximum 
entanglement produced is £o = 9.1 dB and I = 7.5 dB both 
for (3d — ±1.1. Note that as loss in the model is made arbitrar- 
ily small and for a driving field approaching OPO threshold, 
the entanglement becomes arbitrarily strong. 

If we view the total input power to the system as a re- 
source, then it is interesting to see how the strength of en- 
tanglement changes, as the power splitting ratio between the 
fundamental and second-harmonic field is varied, whilst keep- 
ing the total input power held constant. We define the split- 
ting fraction to be i? = |/3i,inP/(5|ai,inP + |/3ianP), and 
follow a path, parameterized by angle On, as defined by 
ad = (/3i,in,c/Q;i,i„,c)y2jsin6'fi, and f3d = V^coaeR. 

In Fig.|4]the total input power ^ has been set to 90% in (a) 
and 400% in (b) of the power required to reach OPO thresh- 
old. A given angle in the plot describes not only the splitting 
fraction for fundamental and second-harmonic power, but also 
the relative phase between them (0 or tt). The radial distance 
in the plot corresponds to the EPR measure of entanglement 
in a dB scale. A circle is drawn at dB to mark the classi- 
cal limit. The shaded areas highlight the entanglement. We 
find that for a total input power required to reach 90% OPO 
threshold, there is a choice of two splittings of i? = 99.9% 
and R = 77%, that both give a maximum amount of entangle- 
ment of £0 — 6.5 dB . These parameters correspond to neither 
a pure OPO, nor a pure SHG process, but correspond rather to 
an OPA operated in a moderately pump-depleted regime. In 
Fig. Ob) we increase the resource available by setting the total 
input power to 400% of OPO threshold. The bi-stability dis- 
played by the system in region (i) is now clearly visible. As 
we drive the system from (i) through to region (ii), the entan- 
glement approaches a strength of Eq = 13.8 dB at i? = 97%. 
The same strength of entanglement is also found at i? = 43% 
which is exactly at the boundary between the regions (iii) and 
(iv). Arbitrarily strong entanglement can be produced in these 



boundary regions, provided that the system has sufficiently 
low loss, and is driven above threshold. 

From these two examples, we can see that stronger en- 
tanglement is made available, when the system is driven by 
higher total input powers. This supports the view, that the to- 
tal input power to the system is a resource for the generation 
of entanglement. To relate this result to a potential experi- 
mental demonstration, we stress that the amount of total input 
power needed to access at least Eq — 6 dB of entanglement is 
near to that required to reach threshold in OPO. Where both 
the detection of such entanglement, and the total input power 
required, are experimentally accessible 0]. 

So far, we have only considered the case where the driv- 
ing fields are coherent states. Although an in-depth discus- 
sion is contained in a forthcoming paper, we can qualitatively 
describe the effect that squeezed driving fields have on the 
harmonic entanglement generation for this system. In short, 
squeezing of the fundamental and/or the second-harmonic 
field in the correct combination of amplitude/phase quadra- 
tures, results not only in stronger entanglement, but also ex- 
tends the regions of strong entanglement closer in toward the 
origin, thereby reducing the total input power needed to ob- 
serve a given strength of entanglement. 

Summary: We have shown that arbitrarily strong entangle- 
ment between a fundamental field and it's second harmonic 
can, in principle, be generated by the second-order non-linear 
interaction. The maximal strength of this harmonic entangle- 
ment, as measured by the EPR and inseparability criteria, is 
produced neither by pure SHG, nor parametric oscillation, but 
rather by an intermediate process. We considered the total in- 
put power that drives the non-linear interaction as a resource 
for the strength of entanglement, and found that an experi- 
mental demonstration of harmonic entanglement using optical 
techniques should be attainable. 
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